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PLANETARY  BOUNDARY  LAYER  PARAMETERIZATION  FOR 
TROPICAL  CYCLONES  BASED  ON  GENERALIZED  SIMILARITY  THEORY 

I.  INTRODUCTION 

The  planetary  boundary  layer  (PBL)  plays  a  critical  role  in 
the  evolution  of  tropical  cyclones  because  the  air-sea  energy  and 
momentum  exchanges  occur  through  the  PBL.  Parameterization  of  the 
PBL  is  thus  a  very  important  aspect  of  tropical  cyclone  modeling. 
Complex  parameterization  is  not  always  computationally  feasible  for 
three-dimensional  or  operational  models.  Computationally  efficient 
yet  accurate  PBL  parameterization  is  imperative  for  these  models. 

The  purpose  of  this  paper  is  to  test  such  a  parameterization  based 
on  the  generalized  similarity  theory  (GST)  in  the  parameterization 
of  the  PBL  in  tropical  cyclone  models. 

The  most  widely  used  PBL  model  is  the  single-layer  parameteri¬ 
zation  where  the  PBL  is  represented  by  one  model  layer  and  the  surface 
fluxes  are  calculated  using  bulk  aerodynamics  formulas.  This  method 
is  simple  and  computationally  economical  but  has  a  major  drawback  in 
that  the  internal  structure  of  the  PBL  cannot  be  resolved.  As  a  con¬ 
sequence,  the  single  layer  PBL  model  using  bulk  aerodynamic  formulas 
can  only  crudely  estimate  various  surface  fluxes. 

An  alternative  approach  where  the  PBL  is  modeled  by  several 
layers  has  been  adopted  in  recent  years  by  Pielke  (1974),  Kurihara 
and  Tuleya  (1974),  Anthes  and  Chang  (1978)  and  Anthes  and  Warner  (1978). 
By  explicitly  resolving  the  PBL  structure,  these  models  often  yield 
more  realistic  estimates  of  the  various  fluxes,  but  the  computational 
cost  is  high. 

Manuscript  submitted  April  1,  1980. 
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Recent  observational  studies  have  facilitated  the  testing  of 
various  PBL  similarity  theories.  The  rather  idealized  assumptions 
characteristic  of  early  theories  have  been  relaxed  and  improved,  and 
more  realistic  and  general  theories  have  been  derived  (Arya  1977; 
Yamada,  1976).  These  developments  have  led  to  a  PBL  formulation 
that  appears  suitable  for  various  PBL  conditions. 

In  the  next  section,  a  brief  review  on  the  generalized  simi¬ 
larity  theory  will  be  given.  In  Section  3,  various  surface  fluxes  in 
typical  tropical  cyclone  conditions  obtained  from  the  GST  will  be  com¬ 
pared  with  those  computed  by  a  time  dependent  one  dimensional  PBL  model 
(Busch,  ef  al . ,  1976;  hereafter  referred  as  BCA) .  The  generalized 
similarity  theory  is  then  applied  to  an  axisymmetric  tropical  cyclone 
model.  The  results  will  be  compared  in  Section  4  with  those  obtained 
using  constant  drag  coefficients. 
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II.  Generalized  Similaritv  Theory 

Similarity  theory  states  that  the  profile  of  any  properl v  scaled 
variable  in  the  PBL  such  as  wind,  temperature,  or  water  vapor,  can  be 
described  by  a  universal  function.  Although  there  has  been  much  debate 
concerning  the  proper  scales  for  the  various  types  of  PBL's,  we  will 
use  the  similarity  theory  that  appears  most  applicable  and  general 
for  tropical  cyclones  (Avra,  1977,  1978;  Yamada,  1976). 

The  mathematical  representation  of  the  PBL  flow  can  be  obtained 
by  assuming  that  there  exists  a  layer  where  both  the  equations  de¬ 
scribing  the  surface  layer  and  the  equations  describing  the  interior 
flow  in  the  mixed  layer  are  both  valid.  In  the  surface  laver,  the 
PBL  variables  are  assumed  to  follow  a  logarithmic-linear  relationship 
with  height 
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where  c  Is  the  unit  vector  parallel  to  the  surface  wind,  P^  Is  the 

turbulent  Prandtl  number,  if  ,  If,  and  1  are  stability  functions,  z 

m  n  q  o 

is  the  roughness  length,  and  0^  and  are  potential  temperature  and 
water  vapor  content  at  z  . 
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For  the  interior  flow  in  the  mixed  layer  we  can  write  the 


equations  in  the  form  of  the  resistance  law 
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where  V  ,  0  and  q  are  the  mean  velocity,  potential  temperatures  and 
mm  m 

water  vapor  content  in  the  PBL,  n  is  a  unit  vector  normal  to  t,  and  h 
is  the  height  of  the  PBL. 

By  matching  (1)  and  (2),  we  obtain  a  set  of  equations  that 
is  simultaneously  valid  for  surface  layer  and  the  mixed  layer. 
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where  u  and  v  are  the  components  of  the  mean  velocity  in  the  PBL  in 
m  m 

t  and  n  directions,  respectively.  Universal  functions  A,  B,  C  and  D 
determined  empirically.  It  is  generally  assumed  that  D  ~  C. 

Using  equation  set  (3),  and  the  Wagara  observational  data, 

Yamada  (1976)  obtained  functions  A,  B  and  C  (eqs.  (13)  through  (18) 
in  his  paper)  which  are  adopted  here  because  of  the  simple  functional 
form  and  accuracy. 

It  should  be  noted  that  V  in  Yamada* s  analysis  is  the  ver- 
tically  averaged  geostrophic  wind  in  the  PBL.  Here  is  the  verti¬ 
cally  averaged  mean  wind.  Thus  the  baroclinic  effects  are  minimized 
(Arya,  1978).  In  addition,  when  actually  applying  Eq.  (3)  in  a  dy¬ 
namic  model,  the  sum  of  the  squares  and  the  quotients  of  Eqs.  (3a) 
and  (3b)  are  used  to  evaluate  the  magnitude  and  direction  of  the 
surface  wind  stress. 

III.  Generalized  Similarity  Theory  and  Multi-Layer  Approaches  in 

a  One  Dimensional  Tropical  Cyclone  Model. 

As  stated  in  the  introduction,  the  Parameterization  of  the  PBL 
based  on  the  GST  for  a  tropical  cyclone  model  is  very  appealing  be¬ 
cause  it  requires  only  one  layer  to  represent  the  PBL.  Thus  it  is 
important  to  test  the  universal  functions  described  by  Eq.  (3)  with 
the  observations  of  tropical  cyclones.  Unfortunately,  detailed  ob¬ 
servations  of  tropical  cyclone  PBLfs  are  not  yet  available,  so  we 
resort  to  use  a  one-dimensional,  multi-layer  PBL  model  for  comparison. 
Table  1  lists  various  models  used  in  this  study. 
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Table  1.  Model  Used  for  Comparison 


Models 

Type 

Dimensions 

PBL  parameterization 

GST 

PBL 

1 

Single-layer,  generalized 
similarity  theory 

BCA 

PBL 

1 

Multi-layer,  Busch  et  al. 
(1976) 

S 

Tropical 

Cyclone 

2 

Single-layer,  generalized 
similarity 

Dl 

Tropical 

Cyclone 

2 

Single-layer,  C  *  0.0015 

cE  =  cH  =  o.oof 

D2 

Tropical 

Cyclone 

2 

Single-laver ,  Cp  =  0.003 

C  E  =CH  =0.003 

The  BCA  1-D  PBL  model  can  reproduce  observations  accurately 
and  it  has  been  used  successfully  in  parameterizing  the  PBL  in  numer¬ 
ical  models  (Anthes  and  Chang*  1978;  Anthes  et  al . ,  1978).  The  model 
has  a  surface  layer  in  which  the  logarithmic-linear  relationships  for 
wind  and  temperature  (Businger  et  al.,  1971)  are  adopted.  It 
should  be  noted  that  the  same  logarithmic-linear  relationship  is 
used  by  Yamada  (1976)  to  generate  the  universal  functions.  Above 
the  surface  layer,  a  time  dependent  mixing  length  is  used.  The 
mixing  length  is  assumed  to  continuously  approach  its  asymptotic 
value,  which  depends  on  the  depth  and  stability  of  the  PBL.  The  e- 
folding  time  of  the  mixing  length  is  proportional  to  the  size  of  the 
energy-containing  eddies  in  the  PBL  and  inversely  proportional  to 
the  turbulent  kinetic  energy. 

The  current  version  of  the  BCA  model  has  20  layers  using  a  stretched 
grid  from  the  surface  to  3  km.  The  lowest  layer  defined  as  the 
surface  layer,  has  a  thickness  of  25  m.  To  properly  simulate 
tropical  cyclone  conditon,  a  rather  strong  pressure  gradient  force 
is  specified  so  that  the  geostrophic  wind  speed  is  strongest  at  the 
surface  and  decreases  with  height  to  90%  of  its  surface  strength  at 
3  km.  The  initial  thermal  state  which  has  been  taken  from  the 
model  hurricane  of  Anthes  and  Chang  (1978)  shows  a  well  mixed  PBL 
below  a  capping  inversion  at  600  m  with  a  sharp  increase  of  potential 
temperature  and  decrease  of  specific  humidity  above  the  inversion. 

The  initial  wind  is  equal  to  the  geostrophic  wind  prescribed  for  each  run. 
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The  BCA  model  was  integrated  12  hours  for  24  cases  with  the 
geostrophic  wind  varying  from  12  to  35  m  s  ^  and  the  sea-surface 
temperature  (SST)  varying  from  24.4  to  29.4°C.  These  variations  pro¬ 
duce  different  stability  conditions  and  large  rates  of  changes  in  the 
initial  hours.  Diagnoses  using  Eq.  (3)  (or  CSA  model)  were  given  at 
model  times  of  1,  2,  3,  4,  6,  8,  10,  and  12  h.  The  vertically  integrated 
u,  v,  0,  and  q  in  the  BCA  model  are  the  mean  PBL  properties  needed 
for  the  diagnoses.  Figure  1  shows  the  comparison  of  the  surface 
wind  stresses  between  the  diagnoses  from  GST  model  and  the  computation 

of  BCA  model.  The  GST  model  produced  nearly  the  same  results  (points 

2 

on  the  diagonal  line  in  Figure  1)  as  the  BCA  model  for  small  u~  value 

and  for  the  unstable  conditions.  The  stress  diagnosed  by  the  GST  are 
higher  than  the  BCA  model  under  stable  conditions.  The  overestimate 
by  the  GST  model  lessens  as  instability  increases. 

The  larger  estimation  of  the  surface  stresses  for  stable  con¬ 
ditions  of  the  GST  model  is  due  to  the  greater  wind  shear  in  the  lower 
PBL  which  is  resolvable  in  the  BCA  model.  This  strong  wind  shear  is 
not  properly  represented  by  the  universal  functions  in  the  single  layer 
approach.  When  such  a  strongly  sheared  layer,  which  usually  accompanies 
a  low-level  inversion,  is  present,  the  PBL  mean  wind  may  not  be  a  proper 
velocity  scale  for  the  PBL.  Fortunately,  the  formation  of  such  a  low- 
level  inversion  is  not  likely  in  tropical  cyclones,  the  GST  model  there¬ 
fore  seems  to  estimate  the  momentum  fluxes  in  typical  tropical  cyclones 
quite  well. 

Figure  2  and  Figure  3  show  the  comparisons  of  the  upward  and 
downward  sensible  heat  fluxes,  respectively.  For  unstable  conditions, 
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the  GST  model  underestimates  the  upward  heat  flux  because  of  the  for¬ 
mation  of  a  strong  superadiabatic  lapse  rate  region  near  the 
surface  in  the  BCA  model.  For  stable  conditions  the  GST  model 
overestimates  the  downward  heat  flux.  This  is  related  to  the  for¬ 
mation  of  low-level  inversion  in  the  BCA  model  which  effectively  cuts 
off  the  downward  transport  of  heat.  These  characteristics  of  the  GST 
model  are  found  in  all  single-layer  parameter izat ions . 

Fig.  4  shows  the  comparison  of  the  moisture  flux  wrq  .  The 

GST  model  gives  a  higher  w  q  for  the  unstable  condition,  and  a 

very  similar  w  q  for  near  neutral  conditions,  and  a 

scattering  in  Fig.  4  is  similar  to  that  in  Fig.  1,  the  more  symmetric 

distribution  about  the  diagonal  in  Fig.  4  is  presumably  related  to 

the  larger  value  of  in  the  lower  part  of  the  BCA  model. 

oz 

These  deviations  of  the  GST  model  from  the  BCA  model  are  much 
less  than  the  model  using  aerodynamic  drag  coefficients.  A  supplementary 
computation  (not  shown)  indicates  various  fluxes  depend  mostly  on  the  wind 
speed  only.  The  deviations  of  the  GST  model  from  the  BCA  model  is 
partially  inherited  from  the  scattering  of  the  original  observation  data 
from  which  the  universal  functions  are  deduced  (Yamada,  1976).  The  compar¬ 
ison,  however,  does  indicate  a  distinct  behavior  of  the  GST  model  for 
different  stability.  It  should  be  pointed  out  the  the  GST  has  been 
applied  to  highly  non-stationary  conditions  in  this  analysis.  We  note 
that  as  the  PBL  approaches  a  quasi-steady  state,  the  GST  model  diagnoses 
becomes  comparable  to  the  BCA  model  results. 
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The  experiments  with  the  one-dimensional  model  indicates 
that  the  GST  model  agrees  very  well  with  the  Busch  model  in  momentum 
and  moisture  flux  predictions.  The  GST  model  does  not  predict  accurate 
fluxes  when  there  are  strong  low-level  gradients  of  potential  temperature. 
The  heat  flux,  however,  is  fortuitously  not  critically  important  in  the 
development  and  maintenance  of  tropical  cyclones  as  pointed  out  by 
Rosenthal  (1971)  and  Anthes  and  Chang  (1978).  The  single  layer  approach 
with  GST  model  in  parameterizing  the  PBL  in  a  tropical  cyclone  is  thus 
adequate. 
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IV.  Incorporation  with  an  Ax i symmetric  Tropical  Cyclone  Model 

a.  Model  Review 

The  axisymmetric  tropical  cyclone  model  used  for  our  inves¬ 
tigation  is  similar  to  the  one  described  in  Anthes  and  Chang  (1978) 
except  that  (1)  the  finest  resolution  at  small  radii  is  30  km,  and 
(2)  there  are  six  vertical  layers.  The  vertical  o-lavers  are  bounded 
at  levels  o  =  0.0,  0.2,  0.3,  0.6,  0.8,  0.93,  and  1. 

Included  in  one  earlier  version  of  the  model  was  a  time 
dependent  equation  for  the  PBL  height  (Deardorff,  1972).  Supplementary 
experiments  showed  that  the  equation  predicts  infinite  growth  of  PBL 
near  the  center  of  the  cyclone  unless  large  horizontal  diffusion  is 
applied.  We,  therefore,  hold  the  height  of  PBu  constant  between  o  = 
0.93  and  1.  The  PBL  has  a  depth  of  650  =  700  m,  a  typical  observed 
PBL  depth  near  topical  disturbances  (Moss  and  Merceret,  1976;  Nitta, 
1974).  Charoock’s  equation  (Delsol  et  al . ,  1970)  is  used  to  compute 
the  roughness  length.  The  initial  conditions  consist  of  a  weak  vortex, 
which  is  in  gradient  balance,  with  maximum  tangential  velocity  of  17 
m  s  \ embedded  in  a  tropical  atmosphere  (Jordan,  1958).  The  lateral 
boundary  conditions  are  Dirichlet  for  the  thermodynamic  variables  and 
zero-vorticity  for  the  momentum. 

For  comparison  purposes  models  were  also  integrated  using 

constant  drag  coefficients  with  C  =  0.0015,  C  =  0.003  (model  Dl) 

D  H 

and  with  C  *  0.003,  C  = 

D  H 


C  0.003  (model  D2). 
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b.  Development  into  Tropical  Cyclones 


The  development  from  the  initial  vortex  into  a 
tropical  cyclone  is  depicted  by  the  maximum  surface  wind  (Fig.  5) 
and  miminum  pressure  (Fig.  6).  After  the  initial  dissipation  due 
to  friction.  Model  S  rapidly  develops  into  an  intense  cyclone  and 
reaches  a  quasi-steady  state  at  about  40  hr.  The  maximum  surface 
wind  is  about  57  m  s  and  the  minimum  surface  pressure  is  approx¬ 
imately  958  mb. 

Models  D1  and  D2  also  intensify  after  the  initial 
dissipation  stage.  These  two  models  deviate  considerably  from 
each  other  and  from  Model  S  after  20  h.  Model  D1  reaches  a  final 
intensity  of  35  m  s  ^and  986  mb,  whereas  Model  D2  reaches  a  final 
intensity  of  43  m  s  ^and  980  mb. 

c.  Cyclone  Structure  at  Quasi-Steady  State 

The  tangential  velocity  field  at  48  h  for  Model  S 
shows  a  strong  and  concentrated  cyclonic  circulation  near  the  center 
(Fig.  7).  The  maximum  velocity  is  at  r  =  30  km  in  the  lowest 
level  and  the  40  m  s  contour  extends  upward  to  o  *  0.6.  The  tangential 
circulation  diminishes  upward  due  to  the  baroclinic  effect  of  the 
warm  core  and  anticvclonic  circulation  occurs  at  the  o  -  0.2  level  out¬ 
side  of  300  km.  The  inflow  is  confined  to  the  PBL  with  maximum  radial 
velocity  of  35  m  s  ^  just  outside  the  region  of  the  maximum  tangential 
velocity  (Fig.  8).  Strong  outflow  with  velocity  more  chan  20  m  s  * 
occurs  at  the  a  m  0.2  level.  Except  for  the  shallow  inflow  and  out¬ 
flow  layers,  the  radial  velocity  is  small  in  mid- troposphere. 


12 


Figure  9  shows  the  relative  humidity  (RH)  at  48  h  of 
Model  S.  The  RH  field  features  (1)  very  moist  PBL  and  outflow 
canopy,  (2)  relative  moist  convective  eyewall  region,  (3)  a  dry  eye 
region  due  to  the  strong  sinking  motion,  and  (4)  a  dry  mid-troposphere 
outside  the  eyewall  due  to  the  general  subsidence. 

The  temperature  field  of  the  quasi-steady  tropical  cyclone 
shows  a  warm  core  with  an  anomaly  of  more  than  +  12°C.  The  +  8°C 
anomaly  extends  to  r  =  300  km  in  the  outflow  layer.  The  +  2^  C 
anomaly  contour  between  a  =  0.6  and  a  =  0.9  outlines  the  sub¬ 
sidence.  The  cooler  PBL  near  the  center  is  due  to  the  adiabatic  expan¬ 
sion  (and  therefore  cooling)  in  the  inflow. 

d.  Air-Sea  Energy  Transfer 

We  will  use  budget  calculations  of  kinetic  energy, 
moisture,  and  heat  to  compare  model  results  and  reveal  the  impor¬ 
tance  of  various  physical  processes.  As  shown  in  Anthes  and  Chang 
(1978),  the  budgets  of  kinetic  energy  (K) ,  moisture  (Q)  and  enthalpy 
(H)  for  a  tropical  cyclone  model  can  be  expressed  schematically  by: 

«  K_  +  C  +  D  +  D  +  Du 
d  v  L  H 

*  Qb  +  E  +  P  +  Qh 

SsHb+A  +  Qp  +  Hs  +  HH 


SK 

dt 


32 

dt 


3H 

3t 
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where  the  terms  with  subscript  b'*-'  denote  tin  lateral  boundary  fluxes, 
and  terms  with  subscript  H's  denote  the  effects  due  to  horizontal 
diffusion.  In  addition,  1)  ,  and  0  represent  the  rates  of  kinetic 
energy  dissipation  due  to  surface  friction  and  cumulus  friction  re¬ 
spectively.  F  and  P  denote  evaporation  and  precipitation  rates  re¬ 
spectively.  0  represents  latent  heat  release,  and  H  the  surface 
♦'-•visible  heat  transfer.  A  denotes  the  adiabatic  heat  in?  rate  due  to 
vertical  moi ion;  (.  denotes  the  kinetic  energy  conversion  rate  due  to 
cr oss-i sobar ic  flow. 

Within  the  model  domain  of  c  -  0  to  1  and  r  s  0  to  300  km, 
the  most  direct  contribution  of  the  PRL  processes  are  the  air-sea  sen¬ 
sible  and  latent  heat  exchanges.  Figure  11  and  12  show  the  time  series 

of  H  and  F.  During  the  early  stages,  the  sensible  and  latent  heat 
s 

transfers  of  D1  are  equal  or  greater  than  those*  of  S.  The  sensible  and 

latent  hea#  exchanges  in  S,  however,  grow  steadilv  after  10  hr  and  finally 

12  -1 

reach  15  x  10  W  and  3  cm  d  respectively,  whereas  in  D1 ,  the  sensible 

1  2 

and  latent  heat  exchanges  only  grow  to  maxima  of  7.5  x  10  W  and  1.8 
cm  d  ^  after  30  h.  After  the  maxima  are  reached,  thev  decrease  with 
time  because  of  stabilization  and  saturation  of  the  PBL. 

Apparently,  the  sensible  and  latent  heat  exchanges  in  S  are  main¬ 
tained  at  high  levels  primarily  due  to  their  nonlinear  dependence  on  wind 
speed  and  stability.  As  wind  speed  in  S  approaches  its  maximum  after  30  h, 
both  and  E  are  enhanced  to  a  greater  extent  than  in  01.  In  addition, 
the  radial  circulation  in  S  is  also  increased  so  that  a  slightly  unstable 
and  dry  PBL  is  maintained  by  the  strong  inflow. 
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The  increased  radial  circulation  is  illustrated  bv  the  con- 


version  rate  (C)  from  the  available  potential  energy  to  kinetic 
energy  shown  in  Fig.  13.  The  conversion  rate  is  approximately  pro¬ 
portional  to  the  combined  strength  of  the  radial  circulation  and  the 
warm  core.  It  is  apparent  that  the  conversion  rate  in  S  grows  rapidly 
between  10  to  20  h,  which  coincides  with  the  increased  evaporation. 

The  proportionally  faster  increase  of  C  than  Hg  and  E  in  Exp.  S  shows 
that  the  increased  heat  transfer  creates  a  stronger  inflow,  which,  in 
turn,  maintains  a  stronger  heat  transfer  and  a  stronger  inward  trans¬ 
port  of  the  dry,  cooler  ambient  air. 

The  two  different  PBL  parameterizations  also  gives  different 
dissipation  rates  (Fig.  14).  Note  that  Dy  in  S  is  smaller  before  15  h 
than  in  Dl ,  and  it  eventually  asymptotes  to  three  times  stronger. 

Again,  this  is  due  to  the  nonlinear  dependence  of  the  surface  stress 
on  wind  speed. 

One  can  argue  that  the  nonlinear  dependence  of  the  air-sea 
energy  exchange  could  be  achieved  by  employing  wind-speed-dependent 
drag  coefficients.  However,  the  parameterization  based  on  similarity 
theory  is  superior.  The  heat  and  momentum  transfers  depend  not  only 
on  wind  speed  but  also  on  PBL  stability,  and  it  is  through  PBL  stability 
that  various  fluxes  are  correlated.  The  surface  roughness  appears  in 
the  formulation  as  a  length  scale  and  thus  the  parameterization 
captures  an  important  feedback  from  the  sea  state.  The  turning  of 
wind  in  low  levels,  neglected  in  the  parameterization  with  drag  co- 
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efficient,  it  also  included  implicitly.  Therefore,  the  surface  stress 


backs  from  the  direction  of  the  PBL  mean  wind,  as  expected. 

Overall,  the  generalized  similarity  theory  appears  to  be 
applicable  in  parameterization  of  the  PBL  in  the  tropical  cyclone 
models.  Such  parameterization  contains  several  important  character¬ 
istics  of  the  PBL  such  as  the  mutual  dependence  of  various  fluxes 
and  the  backing  of  the  low-level  wind,  it  also  produces  very  realistic 
structures  and  energetics  of  tropical  cyclones. 


>i 
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V. 


Summary 


A  PBL  parameterization  based  on  the  generalized  similarity 
theory  has  been  tested  for  tropical  cyclone  models.  We  conclude  that 
a  single  layer  approach  in  parameterizing  the  PBL  is  viable  although 
it  cannot  resolve  the  detailed  structure  of  the  PBL. 

This  parameterization  was  tested  with  the  ID  multi-level  PBL 
model  of  Busch  et  al.  (1976)  under  various  stabilities  and  highly  non¬ 
stationary  conditions.  For  stationary  conditions,  the  computed  water 
vapor  flux,  agrees  well  with  the  Busch  model.  The  momentum  flux  is 
overestimated  and  the  sensible  heat  flux  is  underestimated  especially 
when  there  is  a  low-level  inversion. 

An  axisymmetric  tropical  cyclone  model  which  incorporated 
this  single-layer-parameterization  produced  very  realistic  dynamic 
and  thermodynamic  structures.  In  comparison  with  otherwise  identical 
tropical  cyclone  models  that  employ  constant  drag  coefficients, 
the  superiority  of  stress-dependent  drag  coefficients  in  the  develop¬ 
ment  of  tropical  cyclones  was  demonstrated. 
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Fig,  1  -  Surface  stresses  (ms  )  computed  by  generalized  similarity 
theory  versus  those  computed  by  the  ID  FBL  model  for  stable  (•),  un¬ 
stable  (o),  and  neutral  (m)  conditions.  Points  should  fall  on  the 
diagonal  line  if  two  models  agree. 
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